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ABSTRACT 


\, 

The  interpretation  of  viscous  effects  in  perturbations  of  rapid,  solid- 
body  rotation  of  a  fluid  in  a  cylindrical  gyroscope  is  discussed.  This  is 

s' 

bound  up  with  the  explanation  of  the  roots  of  those  effects  and  aims  to 
improve  the  foundation  on  which  more  sophisticated  questions  in  this  field  can 
be  addressed  successfully.  To  help  with  the  clarification  of  viscous 
displacement  effects,  their  explanation  in  Section  3  is  prefaced  by  a 
discussion  of  two  examples  displaying  some  of  the  salient  features  in  a 
simpler  setting. 
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SIGNIFICANCE  AND  EXPLANATION 


A  number  of  technical  devices,  such  as  some  governors,  navigational 
instruments  and  projectiles,  for  example,  employ  liquid-filled  gyroscopes.  It 
has  long  been  known  that  the  motion  of  the  liquid  can  upset  the  operation  of 
such  a  device  by  means  of  the  forces  and  moments  exerted  by  the  liquid  on  its 
solid  container.  That  fact  has  prompted  studies  of  the  possible  perturbations 
of  the  motion  which  the  fluid  is  meant  to  perform  in  unison  with  its 
container,  first  of  all  under  conditions  in  which  the  contributions  to  the 
forces  and  moments  arising  from  the  viscosity  of  the  fluid  can  be  expected  to 
be  small  in  relative  magnitude. 

A  gyroscope,  however,  is  an  intricate  dynamical  mechanism  that  can  react 
in  rather  startling  ways  to  quite  small  moments,  if  those  are  out  of  step  with 
the  strong,  stabilizing  influence  of  the  spin.  Analysis  of  the  viscous 
effects  shows  that  they  can  act  out  of  step  with  the  main  fluid  motions  and 
therefore  require  more  attention  than  their  magnitude  would  appear  to  justify. 

The  following  notes  respond  to  a  request  for  a  clarification  of  some 
rather  complicated,  conceptual  issues  that  arise  in  the  analysis  of  'small* 
viscous  effects  even  in  the  case  where  the  container  is  a  cylinder  in  steady 
rotation  and  the  fluid  moves  almost  in  unison  with  it.  The  request  reflects  a 
realization  that  a  clearer  understanding  of  the  issues  in  this  simplest  case 
is  needed  as  a  basis  for  fruitful  thought  about  the  still  trickier  issues 
arising  for  different  modes  of  operation  and  for  different  container  shapes  of 
practical  importance  for  such  devices. 

To  divide  the  difficulties,  the  response  is  made  in  three  stages,  of 
which  the  last  addresses  the  viscous  effects  localized  near  the  end-walls  of  a 
spinning  cylinder.  The  first  two  Sections  serve  to  introduce  some  of  the 
salient  issues  in  the  much  simpler  setting  of  two  examples  involving 
oscillation,  but  no  rotation. 
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div  v  -  0,  reduces  exactly  to  3w/3z  »  0  because  nothing  depends  on  x  and  y.  But  on 
the  plate  (envisaged  impermeable),  w  ■  0,  so 

w  i  0  . 

In  the  incompressible  Navier-Stokes  equations,  most  terms  are  now  seen  to  vanish 
identically,  leaving  only 

3u/3t  -  v32u/3z2  ,  (1) 

where  v  is  the  kinematic  viscosity.  The  boundary  conditions  are 


iiut 

u  =  0  e 

at  z  -  0  , 

P 

(2) 

u  ♦  0 

at  z  ♦  •  . 

Initial  conditions  are  not  here  considered,  because  the  simplest  question  of  interest 

concerns  the  motion  of  permanent  character  that  could  develop  ultimately.  The  familiar 
2 

development  variable  z  /(4vt>  is  not  then  relevant;  instead,  the  boundary  conditions  (2) 

show  that  permanent  character  is  possible  only  with  a  motion  periodic  in  time,  and  the 

linearity  of  (1),  (2)  precludes  any  nontrivial  solution  of  frequency  other  than  u. 
i(|)t 

Hence,  u/e  can  depend  only  on  z.  But,  the  problem  has  no  reference  length  other 
1/2 

than  (v/cd)  ,  so  the  motion  can  depend  only  on  the  nondiraensional  distance 

z  (ai/2v>1/2  -  n  (3) 


and  the  velocity  must  have  the  form 


u( z,t)  “  U  f (n )e 
P 


Substitution  of  (4)  into  (1)  yields 


2if  -  f" 


and  the  boundary  conditions  (2)  show  that 

f-e-<1+i>V  (5) 

u  -  «  Re{e'n+i(Ut-n)> 

-n  (6) 

“  U  e  coeOnt  -  n )  . 

P 

This  exact  solution  of  the  Navier-Stokes  equations  is  due  to  Stokes  (1851).  It 
illustrates  two  very  typical  features  of  oscillatory  or  rotating  boundary  layers.  Hie 
first  is  the  appearance  of  a  complex  exponent  in  (5).  Its  real  part  is  not  unexpected  as 
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a  description  of  the  amplitude  decay  of  the  fluid  oscillation  with  increasing  distance 
from  the  plate.  The  imaginary  part  is  seen  from  (6)  to  represent  a  phase  lag  of  the  fluid 
oscillation,  relative  to  the  plate  oscillation,  which  increases  in  proportion  to  the 
distance  from  the  plate.  Such  phase  delay  in  the  fluid  oscillation  turns  out  to  be  a  very 
typical  manifestation  of  the  viscous  diffusion  of  vorticity  in  oscillating  or  rotating 
boundary  layers. 

The  second  feature  typical  of  boundary  layers  caused  primarily  by  oscillation  or 
rotation  is  the  length  scale 

<2v/(i>)1/2  -  « 

characteristic  of  boundary  layer  structure  in  the  direction  normal  to  the  wall.  It  is  the 
basic  boundary  layer  thickness  because  any  specific  choice  of  a  definition  of  such 
thickness  must  clearly  yield  just  a  numerical  multiple  of  6 . 

On  the  other  hand,  Stokes'  very  simple  solution  is  atypical  in  that  here  the  normal 
velocity  w  vanishes  identically  throughout  the  fluid.  Accordingly,  no  definite 
displacement  thickness  can  be  defined,  nor  any  rational  "boundary  layer  thickness"  other 


2.  Second  Example 

For  an  illustration  of  further,  typical  features  in  the  simplest  setting,  envisage 
now  the  same  geometry,  but  with  the  plate  fixed  (and  still  impermeable),  while  a  distant 
contraption  imparts  a  sloshing  motion  to  the  fluid.  It  is  natural  to  expect  such  a  motion 
to  be  non-uniform  in  space,  but  if  the  contraption  is  long  in  one  direction,  the  motion 
should  retain  essentially  the  symmetry  which  makes  the  velocity  independent  of  the  y- 
coordinate  and  makes  the  y-component  of  velocity  vanish  everywhere.  For  present  purposes, 
the  sloshing  motion  can  then  be  specified  adequately  by  saying  that,  if  the  fluid  were 
inviscid,  its  x-velocity  at  the  plate  would  be 

ue(x,t)  «  Us  Uptx/d)  cosmt  (7) 

with  constant,  real  frequency  w,  velocity-amplitude  scale  U8,  and  length  scale  d 
related  to  a  dimension  of  the  distant  contraption  or  the  large  tank  holding  the  fluid. 

For  a  brief  remark,  suppose  the  fluid  were  inviscid,  then  the  mass-conservation 
equation 

3u/3x  +  3w/3z  »  0  (8) 

implies  already  that  w  ?  0,  and  by  (7),  there  must  be  a  vertical  velocity 

we  “  -(UB/d)  uj(x/d)  (*  +  c)  coaut  (9) 

close  to  the  plate.  The  integration  'constant'  c(x,t)  is  zero  if  the  boundary 
condition  w  -  0  at  the  plate  be  imposed  on  this  inviscid  motion,  but  our  interest  in  a 
real  fluid  suggests  a  postponement  of  this  step  until  more  is  learned  about  the  structure 
of  the  transition  to  the  no-slip  boundary  condition. 

Even  for  an  inviscid  fluid,  however,  the  specification  (7)  already  implies  definition 
of  a  non-dimensional  velocity  ratio 

U^/tud)  -  Ro  , 

called  Itossby  number  in  Geophysics.  In  the  present  context,  it  is  natural  to  think  of 
d  as  a  'large'  length  scale,  and  if  the  frequency  id  is  'not  small',  then  the  Rossby 
number  should  have  a  small  value. 
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For  a  Newtonian  fluid,  an  additional,  nondimensional  Ekman  number 

2 

v/(ud  >  -  E 

appears,  which  will  take  even  much  smaller  values  when  the  viscosity  is  'small',  so  that 
the  Reynolds  number 

Re  -  Ogd/v  -  Ro/E 

is  then  large. 

A  rational  way  to  formulate  the  second  example  nondimensionally  is  to  measure 

horizontal  velocity  u  in  units  of  Ug ,  time  in  units  of  u  ' ,  and  horizontal 

distance  x  in  units  of  d,  as  in  (7),  but  vertical  distance  z  again  in  units  of 
1/2 

£  -  (2v/<ii)  to  clarify  the  transition  from  the  sloshing  motion  to  no-slip  on  the 

plate.  Mass-conservation  !8>  then  implies  that  vertical  velocity  w  must  be  measured  in 

units  of  U^j/d,  as  usually  in  boundary  layers.  The  proper  unit  of  pressure,  by 

contrast,  is  not  so  straightforward.  The  first  example  confirmed  the  plausible 

expectation  that  oscillation,  per  se,  does  not  generate  spatial  pressure  gradients,  and  if 

they  occur, they  must  therefore  relate  to  those  of  the  sloshing,  which  can  be  estimated  in 

a  limit  n  ♦  z  ♦  0.  That  limit  can  be  taken  so  as  to  be  governed  plausibly  by  invlscid 

equations,  of  which  the  x-mcmentum  conservation  law  is 

3u  /3t  +  u  3u  /3x  +  w  3u  /3z  -  -p  *3p  /3x  , 
e  e  6  e  e  e 

in  the  dimensional  notation  of  (7)  and  (9).  From  those  two  equations,  u.  +  U  u exp(iut) 

“  I  u 

and  ♦  -u  c( 3u  /3x)exp( iwt)  as  n  *  •,  but  z  ♦  0,  and  accordingly, 

®  s  0 


9pe  iutr  3u0 

3x  "  p“°.«  iu0  *  1,0  3  ( x/d )  1  “0  d  3  ( z/d )  1 


3u 

r  c  0  i  lsiti 

I  un  -  7  TTTaJTI®  f 


irnti 


(10) 


[For  just  a  passing  moment,  u0  is  here  admitted  to  depend  also  on  z.]  For  consistency 
with  the  scaling  of  vertical  velocity  w  just  deduced  from  mass  conservation,  moreover , 
the  constant  c,  if  nonzero,  must  scale  with  5.  Since  present  interest  centers  on  cases 
where  no  and  6/d  are  not  at  all  large,  the  proper  unit  of  pressure  p  is  seen  to  be 
pmdu  • 
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If  all  the  variables  are  now  made  nondiaenslonal  by  reference  to  the  units  just 
specified,  the  mass-conservation  equation  (8)  remains  unchanged  and  the  momentum- 
conservation  equations  of  a  Newtonian  fluid  take  the  form 


3u 

3t 


*  •“(“  fe  ♦  *  If) 


_  -  ifu  i  iin 

3x2  2  3n2 


_  Is.  +  -  ifi  +  1  ifa 

2E  3,1  3x2  2  3n2  ' 


(11) 


(12) 


The  parameters  Ro  and  E  are  seen  to  play  quite  different  roles.  For  'small  viscosity' 
or  high  frequency,  E  ♦  0  and  then  (11)  assumes  boundary  layer  character  and  (12)  shows 
vertical  pressure  variation  across  the  boundary  layer  to  be  negligible,  as  usually,  so 
that  p  ~  p(x,t),  which  is  predictable  from  the  limit  n  ♦  •,  z  ♦  0,  i.e.,  from  the  non- 

dimensional  form  of  ( 10 ) , 


(13) 


where  ue  and  uQ  match  the  notation  in  (7),  (9). 

Comparison  of  (13)  and  (10)  clarifies  why  the  verical  variation  of  the  sloshing 
velocity  plays  no  role  in  the  small-Ekman  limit:  that  variation  is  on  the  length  scale 
d  and  the  analysis  of  the  boundary  layer  is  concerned  with  variations  on  the  length  scale 
5  -  (2E)1//2d,  on  which  the  sloshing  appears  independent  of  x  and  z.  This  explains 
why  the  specification  (7)  of  the  limit  z  ♦  0  of  the  sloshing  velocity  turns  out 
sufficient  for  the  analysis  of  the  boundary-layer  transition.  Furthermore,  the  x- 
dependence  of  the  transition  becomes  a  parametric  effect:  the  transition  depends  only  on 

the  local  values  of  u_  and  3u  /3x. 

•  e 

For  small  Rossby  number  Rs  -  o^/fmd),  moreover,  a  first  approxiaution  to  the 
boundary  layer  structure  is  now  seen  to  be  obtainable  from  a  linear  truncation  of  (11), 
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(14) 


3u/3t  -  U^Su^/St  ■  ”  32u/3n2  , 
almost  as  simple  as  Stokss'  aquation  (1)  for  the  first  example .  me  boundary  conditions 
for  it  appropriate  to  the  description  of  a  transition  from  the  sloshing  (7)  as  n  +  •, 
but  z  ♦  0,  to  no-slip  at  the  plate  n  ■  0  are 
u  ■  0  for  n  ■  0 


u  ♦  u8/Us  -  UQtxJe11  as  1  ♦  “  , 


(15) 


for  all  x,t. 

me  linearity  of  (14),  (IS)  again  rules  out  any  frequency  other  than  u  and  in  fact, 
makes  u/u0  independent  of  x  and  t,  so  that 

U8u(x,n,t)/ue(x,t)  -  e  ifcu(x,n,t)/u0(x)  -  g 
can  depend  only  on  n  and  must  satisfy 

g"  -  2i(g  -  1),  g<0)  -  0,  g(n)  *1  as  n  ♦  •  , 

.  .  -( i+i)n 

whence  g  *■  1  -  e  and 

u  -  uQ(x)eit( 1  -  e-<1+i>n)  .  (16) 

Thus,  [1  -  g(1))u  •  u  -  U  u  exhibits  the  same  diffusive  decay  and  phase  lag  as  f(n) 

®  ®  * 

in  (S). 


me  object  of  the  second  example  is  the  interpretation  of  the  vertical  velocity  w, 
which  is  now  determined  by  (8)  and  the  boundary  condition  w  -  0  at  the  impermeable  plate 
n  ■  0,  so  that 


w(x,n,t) 


n  .  n 

-I  fr  *1*  *  -uo  •lt  /  9(n  •  )dn  > 
*  0 


-»•<*>  eit[n--fT+2-V2e-T’-i<T'+’/4>) 
0  1  ▼  1 


(17) 


An  initially  disturbing  feature  here  is  the  first  term,  which  grows  beyond  bounds  as 
1  ♦  “,  but  actually,  the  first  two  righthand  terms  in  (17)  can  be  recognized  as 
precisely  the  nondimensional  form  of  (9).  me  last  term  describes  an  Internal  boundary 
layer  effect,  which  decays  and  lags  diffusively  with  increasing  distance  from  the  wall. 


and  w  is  seen  to  trail  u  by  another  s/4  in  phase 


The  most  interesting  term  in  the  bracket  of  (17)  is  the  second,  because  it  describes 
a  boundary  layer  effect  which  persists  as  n  +  «  and  even,  for  z  >  0  outside  the 
boundary  layer.  It  describes  a  contribution  to  the  vertical  velocity  which  is  independent 
of  distance  from  the  plate,  but  oscillary  in  time,  and  not  in  phase  with  the  sloshing 
velocity  ufl.  In  the  efficient,  complex  notation  of  (17),  it  is  accordingly  represented 
by  a  complex  multiple  u^(x)/(1  +  1)  of  exp(it),  of  which  the  n-independence  connotes 
persistence  beyond  the  boundary  layer,  the  modulus  describes  amplitude,  and  the  argument 
specifies  phase  difference. 

It  should  be  observed  in  retrospect  that  everything  that  has  1  said  here  about  the 
vertical  velocity  sprang  directly  from  the  mass-conservation  equ.  n  (8).  The  proper, 

physical  interpretation  of  the  effect  under  discussion  is  therefi  in  terms  of  a  mass- 

flow  contribution  to  the  sloshing  motion  arising  from  the  vlscou  .  lion  of  vorticity 
close  to  the  plate.  This  is  the  interpretation  which  goes  to  the  ..cart  of  the  matter. 

On  the  other  hand,  if  one  desires  to  be  done  with  boundary  layer  details  and  to  look 
only  at  their  resultant  effect  in  terms  of  a  corrected  ’inviscid'  description  of  the 
sloshing  motion,  then  the  last  term  in  (17)  becomes  immaterial  and  that  equation  reduces 
to 

iut 

w  -  -(Os/d)u’(x/d)e  [z  -  S/(1  +  i)]  (18) 

in  dimensional,  but  still  complex,  notation.  As  already  noted,  that  is  the  complex 
version  of  (9),  and  there  are  two  alternative,  simple  interpretations.  A  kinematical  one 
is  that  the  plate  z  -  0  appears  porous  to  the  'inviacid  eye'  with  an  oscillatory  suction 
velocity.  That  suction  velocity  is  proportional  to  the  limit  of  dv^/dz  as  z  *  0  and 

it  therefore  shares  the  spatial  structure  of  w#.  its  temporal  structure  is  different, 

-1/2 

however,  because  the  factor  of  -9w^/3z  is  2  S  exp(-ix/4>  >  it  lags  behind  -wp  by 
«/(*»>• 

A  more  common ,  kinematic-geometrical  interpretation  is  in  terms  of  the  displacement 
thickness  4 1  introduced  by  Prandtl ,  probably  in  response  to  early  audiences  sceptical  of 
his  novel  and  controversial  notions  and  particularly,  in  response  to  their  demand  for  a 


concrete  explanation  of  his  order-of-magnitude  concept  of  boundary-layer  thickness.  The 


displacement  thickness  is  defined  as  the  distance  z  ■  4^  from  the  geometrical  plate 
position  at  which  the  corrected  'lnviscid'  description  (18)  of  the  Newtonian  fluid  motion 
satisfies  the  classical  wall-condition  w  «  0.  Hence, 

-  4/(1  +  i)  -  (2v/u)1/2/(1  +  i)  -  (v/u)  1/2e“i1,/4  (19) 

by  (18),  and  this  displacement  thickness  is  complex  because  it  characterizes  a  viscous 
mass-flow  contribution  that  is  out  of  phase  with  the  main  sloshing  motion. 

It  will  be  clear  enough  that  this  is  just  a  simple  example  of  a  mass-flow  effect 

characteristic  of  oscillating  or  rotating  boundary  layers.  Of  course,  it  is  a  small 

1  /2 

contribution  when  4/d  «  (2E)  '  is  a  very  small  number,  but  'small'  must  depend  on  the 

practical  question  at  issue.  Some  issues,  such  as  stability  limits,  are  notoriously 

sensitive  to  phase  relations,  and  the  'small'  viscous  mass-flow  contribution  is  then 

important  because  it  is  the  one  and  only  part  of  the  Newtonian  motion  outside  the  boundary 

layer  that  is  out  of  phase  with  the  imposed  sloshing  (7). 

It  may  be  remarked  finally  that  the  real  part  of  (16)  and  (17)  is 

u  ”  u  Ax, t)  -  a  u  <x/d)e  Z//icos(ut  -  z/4 )  , 
e  SO 


U  1/2  1/2 

cos  **  +  0  cos ( <ut  +  71)  +  (“)  e 


z/4 


coslwt 


1  /2 

4  -  <2v/(ii)  , 

and  illustrates  why  calculations  prefer  to  aim  at  its  simpler,  complex  version  (19).  In 
the  real  version,  the  definition  of  4^  makes  4^  -  -c,  which  fluctuates  between 
positive  and  negative  values. 

That  c  turns  out  independent  of  x  is  seen  in  retrospect  to  have  two  reasons: 
First,  the  analysis  of  the  transition  to  no-slip  concerns  only  the  very  local  dynamics,  in 
the  limit  E  ♦  0,  and  is  influenced  only  parametrically  by  the  x-dependence  of  the 
sloshing,  so  that  ue(x,0,t)  and  -Su^/Sx  at  (x,0,t>  become  the  respective  scale- 
constants  of  the  horizontal  velocity  and  vertical  velocity  gradient  for  the  local 
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dynamics.  Secondly,  the  linearity  of  the  limit  no  ♦  0  makes  the  velocity  components 
directly  proportional  to  their  sealesi  by  suppressing  nonlinear  Interactions,  it 
eliminates  the  more  complicated  structure  of  the  harmonics  of  the  sloshing  frequency  u 
from  consideration. 

It  may  also  be  worth  repetition  that  the  two  mull  parameters  no  and  E  are  seen 
to  play  quite  different  roles  in  the  analysis.  A  double-limit  no  ♦  0,  E  ♦  0  arises, 
and  no  consideration  has  been  given  to  the  manner  in  which  an  experiment  might  approach 
it.  Instead,  the  choice  has  been  made  to  let  no  ♦  0  first  and  then  look  at  (11),  (12) 
for  Ro  -  0  and  E  <<  1.  The  same,  somewhat  arbitrary  choice  is  repeated  in  what 
follows. 
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3.  Ekman  Layer  in  a  Cylinder 

A  question  of  practical  importance  concerns  the  perturbations  of  solid-body  rotation 
of  a  fluid  in  a  spinning  and  nutating  cylinder,  because  of  the  oscillating  forces  and 
moments  which  the  fluid  may  exert  on  its  container.  The  following  focuses  on  the  boundary 
layers  on  the  end-walls  of  the  cylinder,  because  the  obvious  fact  that  they  make  a  very 
small  contribution  to  the  perturbations  at  small  Ekman  numbers  can  be  misleading:  As  in 
the  preceding  examples,  viscous  diffusion  of  momentum  is  associated  with  a  delay,  in  which 
case  the  viscous  contribution  to  the  perturbations  may  be  out  of  phase  with  the  'main' 
fluid  oscillations  and  a  comparison  of  magnitudes  may  be  pointless  in  regard  to  questions 
for  which  phase  relationships  are  decisive.  The  objective  of  the  dicussion  will  be  to 
clarify  why  this  raises  issues  which  are  more  complicated  than,  but  in  their  decieive 
features,  closely  analogous  to,  those  just  discussed  for  the  second  example. 

Envisage,  then,  a  cylinder  of  radius  a  spinning  with  angular  velocity  SI  about  the 

z-axis  and  filled  with  fluid  above  the  plane  end-wall  z  »  0.  For  the  study  of  the 

stability  of  solid  rotation,  it  is  convenient  to  refer  the  motion  to  cylindrical  polar 

coordinates  r,0,z  in  the  frame  of  an  observer  rotating  with  the  cylinder,  who  sees  only 

the  perturbation  from  the  solid  rotation,  but  also,  a  Coriolis  acceleration.  He  can 

absorb  the  centrifugal  force  by  counting  as  his  pressure  p  the  difference  between  the 

1  2  2 

real  pressure  and  the  centrifugal  contribution  -pH  r  +  const  to  it. 

It  is  basic  for  stability  analysis  to  consider  velocity  perturbations  of  scale  Ug 
small  compared  with  the  velocity  fla  of  the  cylinder  (and  such  that  the  space  gradient  of 
velocity  has  only  the  scale  UB/a  and  the  time-rate  of  change  of  velocity,  only  the  scale 

SlU  >.  In  short,  there  is  a  Rossby  number  Ro  •  U  /(aSl>  <<  1  and  to  that  extent,  the 
8  ® 

motion  is  "geostrophlc"  and  a  straightforward  application  of  the  limit  Ro  +  0  will  again 
linearize  the  equations  of  motion  and  will  thereby  ignore  interactions  between 
harmonics.  That  offers  an  opportunity  for  a  relatively  simple,  modal  analysis  by 
separation  of  all  the  variables  so  that  the  resultant  velocity  field  is  the  real  part  of 


i(Ct-B0) 


(20) 


u 

u(r)cos  kz 

e 

V 

m 

v(r)cos  kz 

e 

w 

w(r)sin  kz 

e  J 

where  r,z  are  aade  nondimensional  by  reference  to  a  and  t,  by  reference  to  0  \ 
and  the  velocity  components  are  measured  in  units  of  Ug.  [it  is  assumed,  for  simplicity, 
that  the  height  and  diameter  of  the  cylinder  do  not  differ  so  much  that  their  ratio 
introduces  an  additional  small  parameter.]  The  dependence  on  z  is  left  real  to  display 
that  the  modal  analysis  is  inviscid  to  the  extent  of  recognizing  the  impermeability  of  the 
wall  z  -  0,  but  ignoring  the  no-slip  condition  there. 

In  (20),  m  is  integer  and  C  is  the  nondimensional,  complex  frequency  of  the 
perturbation,  and  it  should  be  remarked  briefly  that  the  analysis  in  the  frame  of  a 
laboratory  -  fixed  observer  differs  only  in  the  interpretation  of  that  frequency.  For 
that  observer,  the  phase  is  constant  when  6  increases  at  an  additional  unit  rate,  so 
that  his  exponent  in  (20)  is  i(C  +  m)t  -  im8  in  the  present  notation. 

The  Newtonian  fluid,  however,  has  also  an  Ekman  number 

E  -  v/(a2fl)  , 

and  when  that  is  small,  all  experience  [Greenspan  1968]  indicates  that  the  transition  to 
no-slip  occurs  again  in  an  Ekman  layer  of  thickness  scale 

«  -  (2v/n>1/2  -  a  (2E)1/2  . 

To  study  that  transition,  requires  again  consideration  of  the  limit 

z  ♦  0,  E  ♦  0  for  fixed  n  -  z/fi  e  (0,»)  .  (21) 

Mass-conservation  shows,  as  in  the  second  example,  that  the  corresponding,  vertical 
velocity  w*  must  be  measured  in  units  iU^/a  -  U8(2E)1//2.  The  issues  raised  by  the 
pressure  scale  are  entirely  analogous  to  those  arising  in  the  second  example:  the  same 
argument  shows  the  scale  of  p  to  be  pfiaU^.  With  these  scale  revisions,  the  Navier- 
Stokes  equations  assume  boundary  layer  character,  and  the  vertical  momentum  equation  shows 
that  pressure  differences  in  the  vertical  direction  within  the  Ekman  layer  can  be  only 
0(E).  flie  limit  (21)  of  the  pressure  perturbation  is  therefore  that  at  z  -  0  in  the 
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inviscid  modal  analysis  loading  to  (20) ,  i.a. 


(22) 


-3p/3r  ■  p[3u#/3t  -  2ve]  +  O(Ro)  , 

z-0 


-r_13p/38  -  p[3v  /3t  +  2u  ]  +  O(Ro) 

0  z-0 


With  horizontal  velocity  components  u,v  still  measured  in  units  of  Ua>  the 
limit  Ro  ♦  0  linearises  the  conservation  equations  for  horizontal  momentum  and  the  limit 
(21)  then  reduces  them  to 


-1  12  2 
3u/3t  -  2v  -  -p  3p/3r  +  —  3  u/3n 

3v/3t  ♦  2u  ■  -p  'r  ^p/Sfl  +  32v/3o2 


(23) 


where  the  pressure  gradient  is  given  by  (22)  and  hence,  depends  only  on  r,8  and  t.  The 

transition  equations  for 

u  -  u  (r,8,0,t>  “  u  and  v  -  v  (r,B,0,t)  •  v 
0  0 

are  therefore  the  typical,  simple  Ekman  equations  of  geostrophic  theory. 


3u/3t  -  2v  -  -j  S^/Sn2 
3v/3t  +  2u  -  a^/Sn2 


(24) 


(25) 


with  boundary  conditions 

u  =  ua  +  u“v5ve  +  v”w’“0  at  n  »  0  , 
u,v  both  +0  as  r\  ♦  •,  but  z  *  0  . 

These  show  again  that  the  Ekman  transition  depends  on  r  and  B  only  paraswtrically 
through  the  effective,  local  scales  u0  and  v0  of  the  linear  system  (24),  (25). 

However,  u  and  v  are  not  directly  proportional  to  ue  and  vfi,  respectively,  because 
of  the  Coriolis  coupling  in  (24).  instead, 

3(u  +  lv)/3t  +  2i(u  +  iv)  »  ~  32(u  +  iv)/3n2 


and  the  boundary  conditions  show 


with 


u  ±  iv  «  (u *  iv >11  -  exp(X  n)J 

0  0  X 


x2  -  2i(C  +  2),  X2  -  21 (C  -  2) 


(26) 
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and  X+,X_ 
unless  both 
here . 


understood  as  the 
Us  C  «  I  2  and 


roots  of  negative  real  part  —  which  are  well-defined, 
la  C  >  0,  and  that  exceptional  case  will  not  be  considered 


In  sun,  the  Oorlolls-coupllnc  Is  left  Intact  by  the  geoatrophlc  Halt  Ho  -  0,  and 
this  gives  a  auch  aore  ccaiplicated  structure  to  the  transition  from  no-slip  to  the 
horizontal  velocity  (ua,v0)  of  the  Bodal  analysis  at  s  -  0.  Bach  of  u  and  v  has 
two  contributions  to  the  transition,  one  of  local  scale  tJ^u^,  the  other,  of  local 
scale  Oava-  Bach  contribution,  aoreover,  displays  a  double  structure,  one  of  boundary 
layer  thickness  4/X+  and  the  other,  of  thickness  4/X_.  The  difference  in  the 
thicknesses,  aoreover,  is  accoapanied  by  a  difference  in  diffusive  phase  lag.  All  the 

saae,  the  transition  shares  with  that  of  the  second  exaaple  the  feature  that  it  is 

1  /2 

coaplete  for  u  and  v,  i.e.,  no  contribution  of  order  B  to  either  is  left  es 

n  ♦  •.  A  mass-flow  effect  of  that  order,  however,  appears  as  possible  here  aa  in  the 

second  exaaple,  and  should  show  up  in  the  noraal  velocity  w. 

Substitution  of  (26)  into  the  mass-conservation  equation 

2  3v  3w 

a7  fru>  ♦  r57  ■  0  (27> 

gives 


-  In  1  n  3w  In  X  n 

n-n  +  i<e  ♦  • "  »  a* <•  ,2e) 

in  the  Unit  (21),  where  stands  for 

1  3 

-757f"e>-73f  '  (29> 

which  is  the  z-coaponent  of  (perturbation)  vorticlty.  the  Ooriolis-coupling  is  seen  frca 
(28)  to  generate  also  two  distinct  contributions  to  w,  one  of  the  scale  Ua%,  just  as 
in  the  second  exaaple,  but  the  other,  of  scale  a  times  the  noraal  consonant  of  the 
perturbation  vorticlty.  The  generation  of  such  vorticlty  by  the  Coriolis  effect  is,  of 
course,  a  well-known  geophysical  phenomenon. 
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Sine*  3w  /3c  and  C  taka  their  values  at  s  »  0  in  the  llait  (21),  tha  first 

a  a 

approximation  to  the  value  of  w  as  h  ♦  ",  but  x  <<  1  is  found  from  (28)  by 
integration  to  be 

-  - 1 .  -  j  mx;1  ♦  x:1)] (»%/».) i „0  ♦  j  isu;1  -  x:1)*.!^  oo) 

and  this  represents  the  part  of  the  normal  velocity  which  persists  beyond  the  boundary 
layer. 


Now,  the  analysis  of  the  Ekman  boundary  layer  is  tailored  to  the  premise  that  the 
subscript  e  distinguishes  quantities  varying  on  the  scales  Q  ,  a  and  Us 
characteristic  of  time,  distance  and  velocity,  respectively,  of  the  core-flow  outside  the 
boundary  layers  on  the  solid  walls  and  outsida  that  region  around  the  rim  of  the  end-walls 
in  which  those  boundary  layers  meet.  Variation  on  those  scales  implies  that  the  (first 
approximation  to  the)  aquations  govsrning  the  core-flow  are  the  limit  Ro  ”  E  »  0  of  (11) 
and  (12).  In  the  cylindrical  coordinates  of  the  rotating  observer,  those  governing 
equations  are  (27)  and 


9ue/3t  -  2v^  -  -3pa/9r  , 

3v^/3t  +  2u^  “  -r  19p#/30  , 

3w^/3t  •  -3p#/3c 

and  the  velocities  (20)  in  the  core-flow  are  effectively  computed  from  them.  Cross- 
differentiation  to  eliminate  the  pressure  from  the  first  two  of  that  equation-trio  yields 

9(rCa>/9t  •  2r9wa/3e  , 

by  (29)  and  (27),  and  since  the  time-factor  of  Ce  is  also  exp(iCt),  by  (29)  and  (20), 

iC  -  <2/C)3w  /3c 

•  ® 


and  (30)  simplifies  to 

w  -  (c  -  j  «(X;1  ♦  X"1)  +  «c“1(X"1  -  X'1))  311^/3*1^0  .  (31) 

This  has  again  the  two  interpretations  noted  in  the  second  example.  The  purely 
kinematical  one  is  that  the  main  effect  of  the  end-wall  transition  on  the  core-flow  is  to 
make  it  look  as  if  there  were  a  small,  very  judicious  amount  of  blowing  and  sucking, 
w  -  w  «  -4(3w  /3c)|  (X  1  +  X  ')  +  C  \X  -  X  ))  , 

m  •  2  ▼  “  "  T 
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through  the  wall  z  -  0,  which  ahares  the  epatlal  structure  of  we,  but  differs  In  phase 
by,  essentially,  the  argument  of  the  complex-valued  square  bracket.  The  more  geometrical 
Interpretation  of  (31)  is  of  a  wall  displaced  to  the  value  z  »  6  at  which  the  normal 
velocity  (31)  vanishes,  which  is 

-  Cj  v/0)1/2((X‘1  +  X"1)  +  (2/OIX"1  -  X"1)]  , 
as  first  determined  by  Kitchens,  Gerber  and  Sedney  [1978].  This  coordinate  z  »  6  1  is 
complex  because  it  is  the  geometric  interpretation  of  a  mass-flow  from  the  end-wall 
boundary  layer  that  is  out  of  phase  with  the  core-flow  structure  of  a  purely  inviscld 
fluid.  That  phase-difference,  in  turn  is  the  reason  why  the  mass-flow  contribution  from 
the  viscous  vorticity-dif fusion  influences  eigenvalues  to  an  extent  out  of  proportion  to 


its  small  magnitude 
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